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Steven J. Schneider 
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Abstract 

The use of annular Hall type MHD generator/accelerator ducts for turbojet energy bypass is evaluated 
assuming weakly ionized flows obtained from pulsed nanosecond discharges. The equations for a 1-D, 
axisymmetric MHD generator/accelerator are derived and numerically integrated to determine the 
generator/accelerator performance characteristics. The concept offers a shockless means of interacting 
with high speed inlet flows and potentially offers variable inlet geometry performance without the 
complexity of moving parts simply by varying the generator loading parameter. The cycle analysis 
conducted iteratively with a spike inlet and turbojet flying at M = 7 at 30 km altitude is estimated to have 
a positive thrust per unit mass flow of 185 N-s/kg. The turbojet allowable combustor temperature is set at 
an aggressive 2200 °K. The annular MHD Hall generator/accelerator is L = 3 m in length with a 
B r = 5 Tesla magnetic field and a conductivity of a = 5 mho/m for the generator and a = 1.0 mho/m for 
the accelerator. The calculated isentropic efficiency for the generator is r| sg = 84 percent at an enthalpy 
extraction ratio, r| Ng = 0.63. The calculated isentropic efficiency for the accelerator is r| sa = 81 percent at 
an enthalpy addition ratio, r| Na = 0.62. An assessment of the ionization fraction necessary to achieve a 
conductivity of a = 1 .0 mho/m is n e /n = 1 .90x 10” 6 , and for a = 5.0 mho/m is n e /n = 9.52x 10“ 6 . 


Nomenclature 


A flow area normal to z axis 

B vector magnetic field 

B 0 tangential magnetic field 

B r radial magnetic field 

B z axial magnetic field 

c e mean random thermal velocity of the electrons 

C p specific heat at constant pressure of air 

E vector of applied electric field 

e magnitude of electron charge 

E 0 tangential applied electric field 

E r radial applied electric field 

E z axial applied electric field 

h enthalpy per unit mass of gas 

I current 

j vector current density 

je tangential current density 

j r radial current density 

j z axial current density 

k Boltzmann’s constant 

K h MHD generator-accelerator Hall loading parameter 

m mass flow 

m e mass of an electron 

MHD Magnetohydrodynamics 

n number density of atoms or molecules in the flow 
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n e electron number density in the flow 

P e Electrical power 

P 0 total pressure of flow 

P s static pressure of flow 

Q electron-molecule momentum transfer cross section 

R gas constant of air 

r radial coordinate of flow 

S entropy 

T e electron temperature 

T 0 total temperature of flow 

U internal energy of gas 

V Hall voltage across generator or accelerator, integral of E z dz 

v vector velocity of flow 

ve tangential velocity of flow 

v r radial velocity of flow 

v z axial velocity of flow 

cox Hall parameter 

z axial coordinate, axis of symmetry of flow 

a fraction of discharge energy going into Joule heating of the flow 

r| Na enthalpy addition ratio of accelerator 

r| Ng enthalpy extraction ratio of generator 

r| sa isentropic efficiency of accelerator 

T| sg isentropic efficiency of generator 

0 tangential coordinate of flow 

p density of gas 

a electrical conductivity 


Introduction 

The Russian AJAX hypersonic vehicle was the first to propose magnetohydrodynamic (MHD) energy 
bypass of a scramjet as a means of extending the scramjet’s performance to higher mach numbers 
(Ref. 1). Analyses of the concept lead to the conclusion that energy bypass of a scramjet can result in 
subsonic ramjet propulsion being maintained in the Mach No. 10 to 16 speed range (Ref. 2). A simplified 
thermodynamic cycle analysis of scramjet energy bypass demonstrated that the concept merits further 
investigation (Ref. 3). Based on these results, an examination of the feasibility of MHD energy bypass 
with turbojets was proposed (Ref. 4). As with the scramjet, the enthalpy into the combustor is to be 
reduced allowing more efficient addition of energy in the combustor without exceeding temperature 
limitations on the turbine materials. Preliminary 1-D analysis of the energy extraction process shows that 
significant enthalpy extraction is possible, but this extraction also results in significant total pressure 
losses (Ref. 5). The annular MHD Hall generator concept was proposed for geometric compatibility with 
the turbojet (Ref. 6). This geometry was successfully applied to in-space propulsion engine by the 
Russians with their Stationary Plasma Thruster using very low pressure xenon plasmas and is flying on 
the Russian Express satellite (Ref. 7). The geometry was also previously explored for stationary 
combustion driven MHD power generation in the 1950s (Ref. 8). A concept vehicle showing power 
generation and distribution is shown in Figure 1. These efforts are motivated by a need to develop a single 
flow path solution to a Mach 7 air breathing engine without the need for mode transitions. 

Besides NASA’s traditional use of plasmas in space thrusters (Refs. 9 and 10), there is a growing 
interest in plasma-based aerodynamics, including flow manipulation through MHD forces, power 
generation, drag reduction, and shock mitigation (Refs. 1 1 to 14). The goal of this paper is to examine the 
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Figure 1 . — Concept vehicle showing MHD energy bypass of a 
turbojet and potential energy management system. 


physics necessary to conduct a 1-D, axisymmetric analysis of an annular MHD generator/accelerator for 
extracting energy electrically by MHD forces from the flow upstream of a turbojet such that the turbojet 
can operate within its thermal limits up to a higher Mach number. This electrical energy can be used for 
vehicle drag reduction or can be reintroduced by MHD acceleration in the nozzle to provide positive 
thrust to the vehicle at the higher Mach number. An important goal is to determine whether 
nonequilibrium ionization efficiently generated in a supersonic flow using pulsed nanosecond discharges 
(Ref. 15) can be used to produce a Lorentz force of a magnitude sufficient to generate substantial amounts 
of electrical power, and to considerably reduce the kinetic energy of the flow. Nonequilibrium ionization 
is believed to be a key factor in this analysis, since the customarily introduced alkali metals are corrosive 
to super alloys. 

The term “plasma” is used to describe a fourth state of matter along with solids, liquids, and gases. 
One usually thinks of them as a high temperature “fully ionized” medium, but the definition of “plasmas” 
is extended to the “partially ionized” gases of interest in this analysis. The most important difference 
between “ionized” and “ordinary” gases is that they can be acted upon by electromagnetic forces. The 
chosen method of analytical treatment called “plasma dynamics” is to define the macroscopic properties 
of the plasma using the theory of the “continuum”. We postulate the macroscopic properties by 
considering the conservation laws of mass, momentum, energy, and charge. In general, however, the 
plasma should be considered as a mixture of ions, electrons, and neutral particles that is neutral overall 
and the “kinetic theory of plasmas” should be applied. The ions and electrons continually experience 
MHD forces which accelerate their motion. All of the particles suffer collisions and the forces acting on 
the ions and electrons are transmitted to the neutrals by collisions thus becoming forces acting on the gas 
as a whole. It is left to future work to investigate the relationship between these two approaches. 

Generalized 3-D MHD Equations 

The fundamental equations of “plasma dynamics” have been greatly simplified by the MHD 
approximation (Ref. 16) which uses the assumptions given below. The assumptions that went into these 
simplified equations (Ref. 17) have been tested experimentally and found to be valid in the Department of 
Energy’s MHD power generation program (Ref. 18) of the 1970s and 1980s. These assumptions are: 

1. The plasma is a single continuous medium of definite composition 

2. Electromagnetic forces are of the same order as gas dynamic forces 

3. The time scale of the problem is the characteristic length divided by a characteristic velocity. 

4. The applied electric field E is of the same order as the induced electromotive force. 

5. The flow velocity is much smaller than the speed of light. 
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6. Maxwell’s equations are unaffected by the gasdynamic motion or, in other words, the magnetic 
field induced by the fluid motion is small compared to the applied magnetic field 

7. Invicid flow is assumed, but a friction term can be easily added to the momentum equation. 

8. No heat loss is assumed, but a heat loss term can be easily added to the energy equation. 

9. The equation of state of the gas is assumed to be the perfect gas law, but other equations of state 
can be implemented. 

10. Body forces due to gravity are neglected 

11. Steady state is assumed 

These simplified 3-D equations (Ref. 17) are as follows: 

Momentum: 

p(v-V)v = j x B - Vp (1) 

Mass: 

V ■ (pv) = 0 (2) 

Energy: 

PvV (^r + u) = —V ■ (vp) + j ■ E (3) 

Current: 

V ■ j = 0 (4) 

Ohm’s Law: 

j = a(E + v x B) — x B (5) 

Axisymmetric 1-D MHD Equations 

The above simplified 3-D equations are reduced to a 1-D axisymmetric form for the annular geometry 
and coordinate system shown in Figure 2 and using the following assumptions. This is similar to a mean 
line flow calculation in a turbomachine. 

1 . Azimuthal symmetry 

2. Large radius approximation, with functions of z only: p(z), v 0 (z), v z (z), j 0 (z), j z (z), T(z), p(z), 

E z (z) 

3. Constant terms: B r , a 

4. Zero terms: B 0 , B z , v r , j r , E 0 , E r 



Figure 2. — Annular geometry 
and coordinate system. 
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Axial Momentum: 



dv z . „ dp 

P V ^ = -)9 Br 'S 

(6) 

Angular Momentum: 

dv 0 . n 

P v z^ = JzB r 

(7) 

Axial Ohm’s Law: 



oB r [E z arcB r ] 

(8) 


Jz - i +(cot) 2 Lb, Vy |b| Vz J 

Angular Ohm’s Law: 

aB r f E z (OTB r 1 

(9) 


Je = .+(«)■ r* “ t ibi + IBI V »J 

Energy: 

? v *dz( V +h )=)z Ez 

(10) 

Continuity: 

pv z A = m 

(11) 

Conservation of Current: 

j z A = I 

(12) 

Electrical power: 

P e = IV 

(13) 

Ideal gas equation of state: 

P = ^ 

K RT S 

(14) 

Calorically perfect gas: 

CpT 0 = C p T s +f + f 

(15) 

Entropy of perfect gas: 

dS=C p — — R— 
P T 0 P 0 

(16) 


1-D Annular MHD Calculation Procedure 

The above axisymmetric 1-D equations describe a pure Hall type device, since the applied tangential 
electrical field E 0 is zero by short circuit around the annulus. In order to numerically evaluate these 
equations, the Hall parameter cox is related to macroscopic properties of the plasma. The following 
definition is used 


GOT = 


e | B | 

m e nQc e 


(17) 


Here we use a two temperature model of the plasma where the electron temperature T e is much higher 
than the bulk gas static temperature T s . Such nonequilibrium plasmas can be generated experimentally in 
high speed flows using high voltage pulsed nanosecond discharges (Ref. 19). Numerical modeling of such 
discharges places the electron temperature at T e = lev, i.e., T e = 1 1605 °K (Ref. 20). The mean random 
thermal velocity of the electrons, c e is calculated to be c e = 6.692xl0 7 cm/sec from 


c e 



(18) 


The electron — molecule momentum transfer cross section for air is estimated (Ref. 17) to be 
Q = 8x1 0“ 16 cm 2 . The number density n for air at the local static pressure P s and local static temperature 
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T s is found from the Loschmidt Number of 2.69xl0 19 molecules/cm 3 at standard temperature and 
pressure. 

The externally applied electric field E z is defined using the Hall loading parameter Kh as follows: 


K h — 


ocE z 

G)TV z B r 


( 19 ) 


where a is the fraction of discharge energy going into Joule heating of the gas. This number has been 
shown to be as low as 0.02 in electrical discharges in nitrogen (Ref. 21). The rest of the discharge energy 
fraction, as high as 0.98, goes into vibrational excitation of the gas and is trapped inside the vibrational 
states of the molecules. As such, it does not contribute to the bulk gas translational temperature which 
also determines entropy rise. The relevant a depends strongly on the gas composition and is especially 
sensitive to the presence of water molecules. Experimental measurements (Ref. 19) indicate that 
a = 0. 1 in dry air and rises to a = 0.4 in moist air. The a in combustion products can be expected to be 
even higher. Energy deposited in vibrational states of molecules is known to have relaxation times on the 
order of milliseconds (Ref. 22), which may be enough time to transit the MHD generator and accelerator 
and affect their performance, but probably not enough time to transit the entire cycle. This effect of a has 
not been considered in the calculations presented in this paper, i.e., the assumed a = 1.0. 

One of the goals of this analysis is to determine the magnitudes of the electrical conductivity, 
magnetic field, and length of flow path that are required for the enthalpy extraction and addition 
operations. These parameters are therefore inputs to the calculation. In order to solve the above set of 
equations, ((6) to (19)), one more parameter such as the flow area must be specified. It was found to 
facilitate the numerical analysis to specify 


d ^ Vz ^ = constant (20) 

instead of specifying the flow geometry. An annular Hall generator/accelerator concept is shown in Figure 3 
where the ionized flow in the annulus interacts with a vector current j composed of components j z and j 0 
which forms a spiral in the annulus. Electrodes providing the pulsed nanosecond discharges for ionization 
are shown placed at several axial positions to maintain the conductivity and the current j z is withdrawn 
through the inlet and exit electrodes and flows through the external power extraction/addition circuit. 



Figure 3. — Annular Hall type MHD power extraction concept for the 
inlet of a turbojet showing the spiral current path through the annulus. 
Power addition uses similar geometry at the exit of the turbojet. 
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MHD Energy Bypass Turbojet Cycle Analysis 

A preliminary cycle analysis was conducted in order to derive gas dynamic inlet parameters for the 
MHD generator and accelerator. It was conducted iteratively with the spike inlet and turbojet until a flow 
path area variation similar to that shown in Figure 4 was attained. A T-s (temperature-entropy) diagram 
for the cycle flying at M = 7 at 30 km altitude is shown in Figure 5. This cycle is estimated to have a 
positive thrust per unit mass flow of 185 N-s/kg. The turbojet allowable combustor temperature is set at 
an aggressive 2200 °K. Note that the inlet total temperature of 2453 °K exceeds this allowable combustor 
temperature. The MHD generator reduces this total temperature to 930 °K by extracting 1.6 MW-s/kg 
from the flow. The turbojet at this altitude is estimated to have only a compression ratio of only 1.09, so 
only 31 kW-s/kg is transferred from the turbine to the compressor. Heat addition in the combustor raises 
the total temperature to 2200 °K and the 1.4 MW-s/kg (generator minus ionizer power) energy bypass to 
the MHD accelerator raises the total temperature of the exit flow to 3633 °K. 

The annular MHD Hall generator is L = 3 m in length with a B r = 5 Tesla magnetic field, a 
conductivity a = 5 mho/m and a loading parameter K h = -0.09. The derivative of the product d(pv z )dz = 
-21.0 kg/m 3 -s is held constant allowing an area variation in the flow path calculation. The Mach number 
reduction from 5.0 to 2.1 and area variation from 0.013 to 0.086 m 2 -s/kg in the generator is shown in 
Figure 6. The total temperature drop from 2530 to 930 °K is shown in Figure 7. This corresponds to an 
enthalpy extraction ratio, r| Ng = 0.63, defined as the ratio of the total enthalpy decrease divided by the 
inlet total enthalpy. This enthalpy extraction results in two orders of magnitude drop in total pressure Po 
from 2350 to 18 kPa as shown in Figure 7. This corresponds to an isentropic efficiency of 84 percent. The 
Hall parameter and axial electric field are given in Figure 8. Keeping the Hall parameter above 20 by 
maintaining low static pressure is important for the enthalpy extraction process indicating that high 
altitude operation is preferable. The total voltage across the 3 m long generator is V = 55.8 kV, i.e., the 
integral of the axial electric field. Total current flow in the generator is I = 28.7 Amp-s/kg. 



Figure 4. — Concept MHD energy bypass engine flow path with spike 
inlet followed by the annular MHD generator, turbojet, and annular 
MHD accelerator. 
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Figure 5. — Total temperature versus entropy for cycle 
flying at M = 7 at 30 km altitude. 
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Figure 6. — Mach number and flow area variation 
in the annular MHD generator. 
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Figure 7. — Total temperature and total pressure drop 
in the annular MHD generator. 
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Figure 8. — Hall parameter and axial electric field 
of annular MHD generator. 


Figure 9. — Mach number and flow area variation 
in the annular MHD accelerator. 


The turbojet is primarily a high altitude combustor since it’s compression ratio is only 1.09. The 
fuel/air ratio needed to raise the temperature to 2200 °K in the combustor is 0.0286 using fuel with a 
heating value of 45000 kJ/kg. Heat addition in frictionless flow in which the flow area is increased in 
such a way that static pressure is held constant in the combustor is employed in the analysis. 

Following a small expansion downstream of the turbine, the annular MHD accelerator receives the 
flow at an inlet Mach Number of M = 1.1, static pressure P s = 58.6 torr and static temperature 
T s = 1790 °K. Based on these conditions, the inlet Hall parameter cox = 51.9. The annular MHD Hall 
accelerator is L = 3 m in length with a magnetic field B r = 5 Tesla, a conductivity a = 1 .0 mho/m, and a 
loading parameter K h = 0.08. The derivative of the product d(pv z )dz = -1.3 kg/m 3 -s is held constant 
allowing an area variation in the flow path calculation. The Mach number increase from 1.1 to 3.1 and 
area variation from 0.069 to 0.094 m 2 -s/kg in the accelerator is shown in Figure 9. The total temperature 
rise from T 0 = 2240 to 3640 °K is shown in Figure 10. This corresponds to an enthalpy addition ratio 
r| Na = 0.62, defined as the ratio of the total enthalpy increase divided by the inlet total enthalpy. This 
enthalpy addition results in a total pressure rise of Po from 17 to 72 kPa as shown in Figure 10. This 
corresponds to an isentropic efficiency of 81 percent. The Hall parameter and axial electric field are given 
in Figure 11. The Hall parameter is above 50 because of the low static pressure and high static 
temperature enabling the efficient enthalpy addition process. The total voltage across the 3 m long 
accelerator is V = 204.9 kV, i.e., the integral of the axial electric field. Total current flow in the 
accelerator is I = 6.7 Amp-s/kg. 

Performance maps of generators and accelerators can be derived by solving Equations (6) through 
(20) for the specified inlet conditions and plotting the exit parameters. When one does this, one finds that 
exit parameters can be plotted against the combination of parameters gLB, indicating that conductivity a, 
length L, and magnetic field B r can be traded off against one another such that the product is constant for 
the same performance. A plot of generator isentropic efficiency r| sg is plotted against gLB for various K h 
in Figure 12 for d(pv z )/dz = -21 kg/m 3 -s. The operating point at r| sg = 0.84 and gLB = 75 on the loading 
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Figure 10. — Total temperature and total pressure rise 
in the annular MHD accelerator. 
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Figure 1 1 . — Hall parameter and axial electric field 
of annular MHD accelerator. 
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Figure 12. — Isentropic efficiency as a function of MHD 
generator loading parameter. 
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Figure 13. — Isentropic efficiency of the MHD generator 
as a function of d(pv z )/dz. 



d(pv z )/dz~kg/m 3 -s 

Figure 14. — Exit flow area as a function of d(pv z )/dz. 
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Figure 15. — Generator power per unit mass flow as a 
function of loading parameter 


curve K h = -0.09 is shown by the circle point. Note that there is no higher performance without increasing 
a or L or B r . A similar plot of generator isentropic efficiency r| sg is plotted against gLB for various 
d(pv z )/dz in Figure 13 for K h = -0.09. Again the operating point is shown by the circle point on the 
d(pv z )/dz = -21 kg/m 3 -s curve. Note that higher performance is available at more negative d(pv z )/dz, but 
this performance comes at the expense of increasing exit area as shown on Figure 14. The power 
generated by the MHD generator is plotted in Figure 15 as a function of gLB and K h for d(pv z )/dz 
=—21 kg/m 3 -s with the operation point shown at 1.6 MW-s/kg. This power minus a fraction used in the 
ionizers is passed on to the MHD accelerator. 
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Figure 1 6.— Accelerator power per unit mass flow as a Figure 1 7.— Isentropic efficiency as a function of MHD 
function of loading parameter K h . accelerator loading parameter. 



Figure 18. — Isentropic efficiency of the MHD accelerator 
as a function of d(pv z )/dz. 


Power used by the accelerator is plotted in Figure 16 as a function of gLB and K h for d(pv z )/dz 
= -1.3 kg/m 3 -s with the operating point of 1.4 MW-s/kg determined by the power passed to it by the 
generator. Operation on the K h = 0.08 loading parameter at gLB = 15 was chosen. A plot of accelerator 
isentropic efficiency r| sa is plotted against aLB for various K h in Figure 17 for d(pv z )/dz = -1.3 kg/m 3 -s. 
The operating point at rj sg = 0.81 and aLB = 15 on the loading curve K h = 0.08 is shown by the circle 
point. Note that higher performance is available at lower K h , but higher a or L or B r is required. A similar 
plot of accelerator isentropic efficiency r| sa is plotted against aLB for various d(pv z )/dz in Figure 18 for 
K h = 0.08. Again the operating point is shown by the circle point on the d(pv z )/dz = -1.3 kg/m 3 -s curve. 
Note that some higher performance is available for less negative d(pv z )/dz, but this component 
performance is limited by an expansion to match the exit static pressure at the altitude of flight. 

An assessment of the electron number density necessary to achieve the required electrical 
conductivity is derived from the following equation: 


a = 


n e e 

m e nQc e 


( 21 ) 


The ionization fraction, n e /n that results in electrical conductivity a = 1.0 mho/m can be determined from 
the equation for a. Substituting Q = 8xl0” 16 cm 2 for the estimated electron-molecule momentum transfer 
cross section for air given previously and c e = 6.692x1 0 7 cm/sec for the mean random thermal velocity of 
the electrons at an electron temperature of T e = lev, i.e., T e =l 1605 °K (Ref. 20), we get 

a = 5.25x10 s — (22) 

n v y 


NAS A/TM— 20 11-217210 


10 


For a conductivity of a = 1 .0 mho/m an n e /n = 1 .90x 10“ 6 is required, and for a = 5.0 mho/m an 

n e /n = 9.52xl0” 6 is required. Ionization fractions were estimated in weakly ionized plasmas using pulsed 

nanosecond, high-voltage discharges sustained by a lower DC voltage (Ref 15). The calculations were 

based on the slope of the DC voltage — current in the sustainer circuit. The reported ionization fractions 

were relatively constant at n e /n = 1.2xl0“ 8 over the static pressure range of 10 to 50 torr of interest here. 

Two to three orders of magnitude higher ionization fraction is required than is currently estimated in the 

measurements. 


Conclusions 

This paper lays a foundation for exploring the thermodynamic benefits of MHD energy bypass in 
hypersonic propulsion with turbojet engines. The equations for a 1-D, axisymmetric MHD 
generator/accelerator which can be used in an integrated combined cycle propulsion system analysis 
incorporating a turbojet are derived and numerically integrated to determine the generator/accelerator 
performance characteristics. MHD power generation and enthalpy extraction in low-temperature, 
nonequilibrium supersonic air flows is explored, with the primary application being an MHD inlet 
diffuser and energy bypass around a turbojet to an MHD accelerator. The annular flowpath geometry of a 
turbojet suggests the desirability of annular type MHD power extraction/addition devices. The concept 
offers a shockless means of interacting with high speed inlet flows and potentially offers variable inlet 
geometry performance without the complexity of moving parts simply by varying the generator loading 
parameter. The 1 -D axisymmetric model of this flow describing the interaction for modeling purposes 
uses the assumptions/approximations of the more generalized 3-D MHD equations for steady-state flow. 
A preliminary cycle analysis was conducted in order to derive gas dynamic inlet parameters for the MHD 
generator and accelerator. It was conducted iteratively with the spike inlet and turbojet flying at M = 7 at 
30 km altitude. This cycle is estimated to have a positive thrust per unit mass flow of 185 N-s/kg. The 
turbojet allowable combustor temperature is set at an aggressive 2200 °K. Performance maps of 
generators and accelerators were derived for the specified inlet conditions. It was found that exit 
parameters can be plotted against the combination of parameters gLB, indicating that conductivity a, 
length L, and magnetic field B can be traded off against one another such that the product is constant for 
the same performance. The annular MHD Hall generator is L = 3 m in length with a B r = 5 Tesla 
magnetic field and a conductivity of a = 5 mho/m. The Mach number reduction is from 5.0 to 2.1 and 
area variation is from 0.013 to 0.086 m 2 -s/kg in the generator. The total temperature drop from 2530 to 
930 °K. This corresponds to an enthalpy extraction ratio, r| Ng = 0.63. This enthalpy extraction results in 
two orders of magnitude drop in total pressure Po from 2350 to 18 kPa. This corresponds to an isentropic 
efficiency of 84 percent. The annular MHD Hall accelerator is L = 3 m in length with a magnetic field 
B r = 5 Tesla and a conductivity a = 1.0 mho/m. The Mach number increase is from 1.1 to 3.1 and area 
variation is from 0.069 to 0.094 m 2 -s/kg in the accelerator. The total temperature rise from T 0 = 2240 to 
3640 °K. This corresponds to an enthalpy addition ratio r| Na = 0 . 62 . This enthalpy addition results in a 
total pressure rise of Po from 17 to 72 kPa. This corresponds to an isentropic efficiency of 81 percent. An 
assessment of the ionization fraction necessary to achieve a conductivity of a = 1.0 mho/m is n e /n 
= 1.90xl0“ 6 , and for a = 5.0 mho/m is n e /n = 9.52x1 0" 6 . 

Future Work 

Cycle analysis using this approach should continue with a computational method set up to rapidly 
establish the performance maps of the generator and accelerator as a function of inlet conditions and 
control parameters. The method should iterate with the performance of a turbojet to establish and 
optimize the cycle performance as a function of flight regime. The magnetic field B r , the load factor K h , 
and the conductivity a are all variables that can be varied without any change in geometry. Mass bypass 
around the turbojet may prove beneficial and can be incorporated with geometric change. The analysis 
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should follow the established dynamic pressure corridor of the airframe. Once a 1-D axisymmetric design 
of a generator and accelerator is established, the large radius approximation should be relaxed going into 
a 2-D axisymmetric analysis to calculate hub to tip variations in the flow. Finally a full 3-D analysis 
should be conducted, allowing for 3-D variations in the vectors B, E, a, and j. 

Since the conductivity of the working fluid in the cycle analysis is based on a two temperature gas 
model with a nonequilibrium ionization process, a higher fidelity analysis should be established in which 
the plasma is considered as a mixture of ions, electrons, and neutral particles that is neutral overall. This 
analysis would use the “kinetic theory of plasmas” to follow the relaxation time of the high temperature 
electrons and the resulting loss of conductivity. The analysis should also look at energy transfer from the 
electrons into the vibrational states of the neutrals and their subsequent relaxation times. 

Experimentally, the pulsed nanosecond discharge method of ionization needs to be developed with 
the goal of delivering a higher ionization fraction than that currently estimated. Methods of measuring this 
ionization fraction and the attendant electron temperature distribution function need to be established. 
Also, methods of measuring the fraction of discharge energy a going into Joule heating of the gas need to 
be established along with the variability of a as a function of gas composition, especially air with water 
vapor and carbon dioxide present. Proof of concept generator and accelerator tests need to be designed for 
wind tunnel testing. These tests should be conducted in parallel with the analysis which would seek to 
reproduce experimental results analytically. 
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